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Abstract 

Recent advances in string theory have highlighted the need for reliable numerical meth- 
ods to calculate correlators at strong coupling in supersymmetric theories. We present a 
calculation of the correlator (0|T ++ (r)T ++ (0)|0) in N = 1 SYM theory in 2+1 dimen- 
sions. The numerical method we use is supersymmetric discrete light-cone quantization 
(SDLCQ), which preserves the supersymmetry at every order of the approximation and 
treats fermions and bosons on the same footing. This calculation is done at large N c . 
For small and intermediate r the correlator converges rapidly for all couplings. At 
small r the correlator behaves like 1/r 6 , as expected from conformal field theory. At 
large r the correlator is dominated by the BPS states of the theory. There is, however, 
a critical value of the coupling where the large-r correlator goes to zero, suggesting 
that the large-r correlator can only be trusted to some finite coupling which depends 
on the transverse resolution. We find that this critical coupling grows linearly with the 
square root of the transverse momentum resolution. 
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1 Introduction 



Our original motivation to study correlators of the energy momentum tensor |l], was 
the discovery that certain field theories admit concrete realizations as a string theory on 
a particular background ||. Attempts to apply these correspondences to study the details 
of these theories have only met with limited success so far. The problem stems from the 
fact that this correspondence relates weakly coupled supergravity and strongly coupled SYM 
theory. Unfortunately we only have firm control of either theory in the weak coupling limit. 
The objective of our program is to improve this situation substantially. 

Previously we showed that Supersymmetric Discrete Light Cone Quantization (SDLCQ) 
[|. BJ can be used to solve supersymmetric field theories in the strong coupling limit [[], |7|, |J. 
This then allowed us to make a quantitative comparison between the strongly coupled N = 
(8, 8) SYM theory and the supergravity approximation of the string theory |l], |2j in 1+1 
dimensions. The SDLCQ approach works particularly well in 1+1 dimensions; however, 
it can be extended to more dimensions. Recently, we solved for the spectrum and wave 
functions of M = 1 SYM in 2+1 dimensions || llOfl . 

Aside from our numerical solutions, there has been very little work on solving SYM the- 
ories using methods that might be described as being from first principles. While selected 
properties of these theories have been investigated, one needs the complete solution of the 
theory to calculate the correlators. By a "complete solution" we mean the spectrum and 
the wave functions of the theory in some well-defined basis. The SYM theories that are 
needed for the correspondence with supergravity and string theory have typically a high 
degree of supersymmetry and therefore a large number of fields. The number of fields sig- 
nificantly increases the size of the numerical problem, and, therefore, in this first calculation 
of correlators in 2+1 dimensions we consider only M = 1 SYM. 

A convenient quantity that can be computed on both sides of the correspondence is 



the correlation function of a gauge invariant operator [11, O]. We will focus on two-point 



functions of the stress-energy tensor. This turns out to be a very convenient quantity to 
compute for reasons that are discussed in H. Following the procedure that we used in 
our calculation in 1+1 dimensions |l|, 0, we continue the results to Euclidean space. The 
correlator of the energy momentum operator has been studied in conformal field theory in 



2+1 dimensions [|13[], and this provides a reference point for our results. The structure of 
the correlators in conformal field theory is particularly simple in the collinear limit x± — > 0, 
and we therefore find it convenient to work in this limit. From results in conformal field 
theory we expect that correlators behave as 1/r 6 at small r, where we are probing deep 
inside the bound states. We have confirmed this 1/r 6 behavior by an analytic calculation of 
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the free-particle correlator in the DLCQ formalism |I| . 



The contributions of individual bound states have a characteristic length scale corre- 
sponding to the size of the bound states. On dimensional grounds one can show that the 
power behavior of the correlators are reduced by one power of r; so for individual bound 
states the correlator behaves like 1/r 5 for small r. It then becomes a nontrivial check to 
see that at small r the contributions of the bound states add up to give the expected 1/r 6 
behavior. We find this expected result as well as the characteristic rapid convergence of 
SDLCQ at both small and intermediate values of r. 

At large r the correlator is controlled by the massless states of the theory. In this theory 
there are two types of massless states. At zero coupling all the states of the 1 + 1 dimensional 
theory are massless, and for non- vanishing coupling the massless states of the 1 + 1 theory are 
promoted to massless states of the 2 + 1 dimensional theory fllQf . These states are BPS states 



and are exactly annihilated by one of the supercharges. This is perhaps the most interesting 
part of this calculation because the BPS masses are protected by the exact supersymmetry 
of the numerical approximation and remain exactly zero at all couplings. Commonly in 
modern field theory one uses the BPS states to extrapolate from weak coupling to strong 
coupling. While the masses of BPS states remain constant as functions of the coupling, 
their wave functions certainly do not. The calculation of the correlator at large r provides 
a window to the coupling dependence of the BPS wave functions. We find, however, that 
there is a critical coupling where the correlator goes to zero, which depends on the transverse 
resolution. A detailed study of this critical coupling shows that it goes to infinity linearly 
with the square root of the transverse resolution. Below the critical coupling the correlator 
converges rapidly at large r. One possible explanation is that this singular behavior signals 
the breakdown of the SDLCQ calculation for the BPS wave function at couplings larger 
than the critical coupling. If this is correct, calculation of the BPS wave function at stronger 
couplings requires higher transverse resolutions. We note that above the critical coupling 
(see Fig. ||] below) we do find convergence of the correlator at large r but at a significantly 
slower rate. 

This paper is organized as follows. In section 2 we discuss light cone quantization and 
SDLCQ. The correlators are discussed in section 3 for the free theory and in section 4 for 
the full theory. In section 5 we discuss our numerical results. A brief conclusion is given in 
section 6. 



2 



2 Light-Cone Quantization and SDLCQ 



The technique of DLCQ is reviewed in [14], so we will be brief here. The basic idea of light- 
cone quantization is to parameterize space-time using light-cone coordinates 
and to quantize the theory such that x + plays the role of a time. In the discrete light-cone 
approach, we require the momentum p_ = p + along the x~ direction to take on discrete 
values in units of P + /K, where P + is the conserved total momentum of the system and 
K is an integer usually referred to as the harmonic resolution |TJ]]. One can think of this 
discretization as a consequence of compactifying the x~ coordinate on a circle with a period 
2L = 2nK/P + . Along the direction x L the transverse momentum is discretized as well; 
however, it is treated in a fundamentally different way. The transverse resolution is T, and 
we think of the theory as being compactified on a transverse circle of length I. Therefore, 
the transverse momentum is cut off at ±27rT/l and discretized in units of 2ir/l. Removal of 
this transverse momentum cutoff therefore corresponds to taking the transverse resolution 
T to infinity. 

The advantage of discretizing on the light cone is the fact that the dimension of the 
Hilbert space becomes finite. Therefore, the Hamiltonian is a finite dimensional matrix, and 
its dynamics can be solved explicitly. In SDLCQ one makes the DLCQ approximation to the 
supercharges, and these discrete representations satisfy the supersymmetry algebra. There- 
fore SDLCQ enjoys the improved renormalization properties of supersymmetric theories. Of 
course, to recover the continuum result we must send K and T to infinity and, as luck would 
have it, we find that SDLCQ usually converges faster than ordinary DLCQ. Faster conver- 
gence is important because the size of the matrices and, consequently, the difficulty of the 
computation grow as the resolution is increased. 

Let us now review these ideas in the context of a specific super- Yang-Mills theory. We 
start with 2 + 1 dimensional Af = 1 super- Yang-Mills theory defined on a space-time with 
one transverse dimension compactified on a circle. The action is 

r l , ,1 



S = Jd 2 xJ Q dx ± ti(--F'" / F fU/ + itt 7 "D M tt). (L 



After introducing the light-cone coordinates x ± = ^75 i^ 1 ), decomposing the spinor \1> 
in terms of chiral projections 

1 + 7 5 1 - 'Y 5 

* = -2Vr*> x = ( 2 ) 

and choosing the light-cone gauge A + = 0, we obtain the action in the form 

"I 



S = dx dx / dx 1 tr 



(d^A-y + D +( f)d^ + iipD+iP+ 
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(3) 



A simplification of the light-cone gauge is that the non-dynamical fields A and x may be 
explicitly solved from their Euler-Lagrange equations of motion 



A' 



1 



d 2 



(P. 



V2d^ 



D±ip. 



(4) 



These expressions may be used to express any operator in terms of the physical de- 
grees of freedom only. In particular, the light-cone energy, P~ , and momentum operators, 
P + ,P X , corresponding to translation invariance in each of the coordinates x and x± may 
be calculated explicitly as 



P" 



dx / dx±tr 
Jo 

dx~ / dx±tr 
Jo 

i 



P± = J dx J dx±tr [d-(pd±(f) + iipd±ip] . 



(5) 
(6) 
(7) 



The light-cone supercharge in this theory is a two-component Majorana spinor, and may be 
conveniently decomposed in terms of its chiral projections 



Q + = 2 1/4 J dx~ J* dx ± ti [<j>d-ij) - ipd. 
Q- = 2^Jdx~J 



dx i tr 



2c\<^ + 5 ym (i[0, <M] + 2^) — $ 



(8) 
(9) 



The action gives the following canonical (anti-)commutation relations for propagating 
fields for large N c at equal x + : 

<frij(x~,x ± ),d-<p k i(y~,y ± )] = {ipij(x~ , x±),ip k i(y~ ,V±)} = \ 8 { x ~ ~ V~) s (%± ~ V±)$u&jh- 

(10) 

Using these relations one can check the supersymmetry algebra 



{Q+, Q+} = 2V2P+, {Q-,Q-} = 2V2P-, {Q + , Q-} = -4Pi 



In solving for mass eigenstates, we will consider only states which have vanishing trans- 
verse momentum, which is possible since the total transverse momentum operator is kinemat- 
ical.[] On such states, the light-cone supercharges Q + and Q~ anti-commute with each other, 

1 Strictly speaking, on a transverse cylinder, there are separate sectors with total transverse momenta 
2ttN±/L; we consider only one of them, N± = 0. 
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and the supersymmetry algebra is equivalent to the Af = (1,1) supersymmetry of the di- 
mensionally reduced (i.e., two-dimensional) theory Moreover, in the P± = sector, the 
mass squared operator M 2 is given by M 2 = 2P + P~ . 

As we mentioned earlier, in order to render the bound-state equations numerically tract- 
able, the transverse momenta of partons must be truncated. First, we introduce the Fourier 
expansion for the fields and ip, where the transverse space-time coordinate x± is periodically 
identified 



4>ij(0,x ,x± 

1 



E 



dk+ 



ipij(0,x~,x ± ) = 



CLij(k ,n )e 



±\-ik+x'-i^^-x ± , t (i + 



— X ^ +a) i (k + ,n ± )e 



E 



2V7r/ ± — Jo 



dk~ 



k+x~ 



-1 ; X± 



+ bUk + 1 n ± )e [k+x ' +[2 -^^ 



(13) 



Q- 



+ 



Substituting these into the (anti-)commutators fllQl) , one finds 

a ij (p + ,n_ L ),a] k (q + ,m_ L )'\ = {%(p + , n±), b\ k {q + , m±)} = S(p + - q + )5 n±im± 5u5 jk . (12) 
The supercharges then take the following form: 

POO r 

Q+ = 12 1 / 4 J2 / dkVk \b\jik, n^aijik, n x ) - a^k, n^b^k, n x ) 

i E / dk—r= aUk,n )bij(k,n ) - b\Ak,n )a ij (k,n ) 

E J Q dk 1 dk 2 dk 3 6(k 1 + k 2 - k 3 )S n ± +n ± in ±S^ 

[ a lk( k i, n^)alj(k 2 , n^bijih, 713 ) - blj(k 3 , 713)0^(^1, n±)a kj (k 2 , n£)] 
[al k (k 3 ,nj)a kj (k 1 ,n^)b ij (k 2 ,n2) - a} fc (fci, n^&J^Afe, n 2 )aij(k 3 , )] 
■[^(^1)^1)4^(^2,^2)0^(^3,^) - a\j{k 3: n^)b ik (ki)a kj (k 2 , n^)] 
: (h, n^)blj(k 2 , rig )&y(A&, rig ) + b\j(k 3 , n^)b ik {k x , n x )b kj (k 2 , n£)] 



£-1/4 



9ym 



n 1 £Z 
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k 2 
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2VA;iA; 2 
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2^3 
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+ k 3 


2Vhh 




ki 



1 1 



(14) 



We now perform the truncation procedure; namely, in all sums over the transverse momenta 
nf appearing in the above expressions for the supercharges, we restrict summation to the 
following allowed momentum modes: nj- = 0, ±1... ± T. Note that this prescription is 
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symmetric, in the sense that there are as many positive modes as there are negative ones. In 
this way we retain parity symmetry in the transverse direction. The longitudinal momenta 
hi = riiir/L are restricted by the longitudinal resolution according to K = J2i n i- 

There are three commuting Z2 symmetries. One of them is the parity in the transverse 
direction, 

P : aijikjU 1 ) -> -aij(k, -n 1 ), bi,j(k,n L ) -> b^k, -n 1 ). (15) 
The second symmetry [|l5j is with respect to the operation 

S : a i j(k,n ± ) — > —aj^k,^), b i j(k,n ± ) — > —bji(k,n L ). (16) 

Since P and S commute with each other, we need only one additional symmetry R = PS 
to close the group. Since Q~ , P and S commute with each other, we can diagonalize them 
simultaneously. This allows us to diagonalize the supercharge separately in the sectors with 
fixed P and S parities and thus reduce the size of matrices. Doing this one finds that the roles 
of P and S are different. While all the eigenvalues are usually broken into non-overlapping 



^-odd and 5*-even sectors |16|], the P symmetry leads to a double degeneracy of massive 



states (in addition to the usual boson- fermion degeneracy due to supersymmetry) . 



3 Free particle Correlation Functions 



Let us now return to the problem at hand. We would like to compute a general expression 
of the form 

F(x + ,x-,x ± ) = (0|T ++ (x + ,x-,x ± )T ++ (0,0,0)|0). (17) 

Here we will calculate the correlator in the collinear limit, that is, where x x = 0. We know 
from conformal field theory [HJ calculations that this will produce a much simpler structure. 

The calculation is done by inserting a complete set of intermediate states \a), 

F{x + ,x~,x ± = 0) = ^(0|T ++ (x-,0,x ± = 0)|a)e- iP " :c+ (a|r ++ (0,0,0)|0). (18) 



with energy eigenvalues P a . In |T7| we found that the momentum operator T ++ (x) is given 

by 



T ++ (x) = tr 



(<9-0) 2 + - (iipd.ip - i{d_i))i)) 



T+ + (x)+T+ + (x). 



(19) 



In terms of the mode operators, we find 



T ++ (a; + ,x-,0)|0) = — £ £ T(n,m)e 

£± J '> n,m n±,mx 



-i{P++P+)x~ 



|0> 



(20) 
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where the boson and fermion contributions are given by 

L I . . . yj Til 72 r x \t/ 

— TgMn, m) 0) = — - — tr a-,-(n, njja L(m, mjj 0) 
71 2 L J 

and 



L^,i , , , , (n — m) 
-T£ + {n,m)\0) = - — ^tr 

7T ' 4 



6i,-(n,njL)65 i (m,m±) |0) . 



(21) 



(22) 



Given each \a), the matrix elements in (18) can then be evaluated, and the sum computed. 

First, however, it is instructive to do the calculation where the states |a) are a set of free 
particles with mass m. The boson contribution is 

( it \ 2 

n,m,s,t 

X x'llllis/ ( ' -",„■■ 



F(x + ,x-,0) B = WfTi) (0\tT[a(n } n ± )a(m,m ± )}ti[a\s,s ± )a\t,t ± )}\0} (23) 



where the sum over n implies sums over both n and n±, and 

m 2 + (2n ± vr/Z) 2 



p: 



and P„ 



nir 



2mr/L ~ n L 

The sums can be converted to integrals which can be explicitly evaluated, and we find 



(24) 



F{x + ,x-,0) B 



: m b — 



2(2vr^ 

where x 2 = 2x~x + . Similarly for the fermions we find 

2 



X X 



K 2 /2 (mx) 



(25) 



F(x + ,x- J 0) F = Y, (vkl) (0|tr[6(n,n±)6(m,m±)]tr[6t(a,s ± )6t(t,tj.)]|0> (26) 



t \8LHj 
x (m — n)(s — t) e 

After doing the integrals we obtain 



F(x + ,x",0) F 



i 5 ( x + \ 1 



4(2tt) 



-m 



x x 



K 7 / 2 (mx)Kn/>>.(mx) - K K/0 (mx 



-3/2 



5/2V 



(27) 



We can continue to Euclidean space by taking r = \/2x + x~ to be real, and, finally, in the 
small-r limit we find 



x+) F ^ +,X ' 0) 8(2vr) 2 r6 



(28) 



which exhibits the expected 1/r 6 behavior. 
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4 Correlation Function in SDLCQ 



Now let us return to the calculation using the bound-state solution obtained from SDLCQ. 
It is convenient to write 

F(x\x~,0) = £ (^) 2 (0|-T(n,m)e-^> + -^-T( S ,t)|0), (29) 

where P~ p is the Hamiltonian operator. We again insert a complete set of bound states \a) 
with light-cone energies p- = (Ml + Pl)/P + at resolution K (and therefore P + = nK/L) 
and with total transverse momentum P± = 2Nj_n/l. We also define 

\u) = N u — V 5 n+rn , K 5n ± +m ± ,N ± T(n, m) |0) , (30) 

7T — 

" n,m 

where N u is a normalization factor such that (u\u) — 1. It is straightforward to calculate 
the normalization, and we find 

1 K 3 , 1 w 
_ = _ (1 --)(2T + 1). (31) 

The correlator (E9l) becomes 



F(X+ ' X "'° )= ^ (2T^) 2e " iP " + " lP+ "iv^ IHa)|2 - (32) 

We will calculate the matrix element (u\a) at fixed longitudinal resolution K and trans- 
verse momentum iVi = 0. Because of transverse boost invariance the matrix element does 
not contain any explicit dependence on Nj_. To leading order in 1/K the explicit dependence 
of the matrix element on K is K 3 ; it also contains a factor of /, the transverse length scale 
To separate these dependencies, we write F as 

F(x + x~0)-- V 1 1 ( 7rK ) 3 c -iP^-iP+*- 1W1 2 f33 } 

We can now do the sums over K and N± as integrals over the longitudinal and transverse 
momentum components P + = ttK/L and P L = 2nN±/l. We obtain 

(0F^, X ~, O) ^^M^K M ^ (34) 



In practice, the full sum over a is approximated by a Lanczos |L8| iteration technique 
that eliminates the need for full diagonalization of the Hamiltonian matrix. For the present 
case, the number of iterations required was on the order of 1000. 



Looking back at the calculation for the free particle, we see that there are two independent 
sums over transverse momentum, after the contractions are performed. One would expect 
that the transverse dimension is controlled by the dimensional scale of the bound state 
(Rb) and therefore the correlation should scale like l/r 4 R%. However, because of transverse 
boost invariance, the matrix element must be independent of the difference of the transverse 
momenta and therefore must scale as l/r 5 i?^. 



5 Numerical Results 

The first important numerical test is the small-r behavior of the correlator. Physically we 
expect that at small r the bound states should behave as free particles, and therefore the 
correlator should have the behavior of the free particle correlator which goes like 1/r 6 . We 
see in ([33]) that the contributions of each of the bound states behaves like 1/r 5 . Therefore, 
to get the 1/r 6 behavior of the free theory, the bound states must work in concert at small r. 
It is clear that this cannot work all the way down to r = in the numerical calculation. At 
very small r the most massive state allowed by the numerical approximation will dominate, 
and the correlator must behave like 1/r 5 . To see what happens at slightly larger r it is useful 
to consider the behavior at small coupling. There, the larger masses go like 

M^T. { M- (35) 

i 

Consequently, as we remove the k x cutoff, i.e. increase the transverse resolution T, more and 
more massive bound states will contribute, and the dominant one will take over at smaller 
and smaller r leading to the expected 1/r 6 . This is exactly what we see happening in Fig. [I] 
at weak coupling with longitudinal resolution K = 4 and 5. The correlator converges from 
below at small r with increasing T, and in the region —0.5 < logr < 0.5 the plot of r 5 
times the correlator falls like 1/r. In Fig. |2| at resolution K = 5 we see the same behavior 
for strong coupling (g = g YM ^N c l/27r 3 ^ 2 = 10) but now at smaller r (logr ~ —0.5) as one 
would expect. Again at strong coupling we see that the correlator converges quickly and 
from below in T. All indications are that at small r the correlators are well approximated by 
SDLCQ, converge rapidly, and show the behavior that one would expect on general physical 
grounds. This gives us confidence to go on to investigate the behavior at large r. 



The behavior for large r is governed by the massless states. From earlier work ||, |10 



on 



the spectrum of this theory we know that there are two types of massless states. At g = 
the massless states are a reflection of all the states of the dimensionally reduced theory in 
1 + 1. In 2+1 dimensions these states behave as g 2 M 2 +1 . We expect therefore that for g — 
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3 



3 




vs. r (a) in units where g = (7ymv^VcI/27t 3//2 = 0.10 for K = 4 and T = 1 to 9; (b) in units 



where g = ^ Y mv / ^/2vt 3/2 = 1 for K = 5 and T = 1 to 9. 
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-6 H i i i i i i i i i i i i i i i ' ' ' ' i ' ' ' i— 
-3-2-1012 

logio r 



Figure 2: Same as Fig. |(b), but for g = ^ Y mvW/2tt 3 / 2 = 10. 
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K = 5 and (b) K = 6 and for various 



there should be no dependence of the correlator on the transverse momentum cutoff T at 
large r. In Fig. |l|(a) this behavior is clearly evident. 

At all couplings there are exactly massless states which are the BPS states of this theory, 
which has zero central charge. These states are destroyed by one supercharge, Q~ , and not 
the other, Q + . From earlier work || on the spectrum we saw that the number of BPS states 
is independent of the transverse resolution and equal to 2K — 1. Since these states are exactly 
massless at all resolutions, transverse and longitudinal convergence of these states cannot 
be investigated using the spectrum. These states do have a complicated dependence on the 
coupling g through their wave function, however. This is a feature so far not encountered 
in DLCQ |14|. In previous DLCQ calculations one always looked to the convergence of the 
spectrum as a measure of the convergence of the numerical calculation. Here we see that 
it is the correlator at large r that provides a window to study the convergence of the wave 
functions of the BPS states. In Fig. |2| we see that the correlator converges from above at 
large r as we increase T. 

We also note that the correlator at large r is significantly smaller than at small r, par- 



ticularly at strong coupling. In our initial study of the BPS states [[l(J we found that at 
strong coupling the average number of particles in these BPS states is large. Therefore the 
two particle components, which are the only components the T ++ correlator sees, are small. 

The coupling dependence of the large-r limit of the correlator is much more interesting 
than we would have expected based on our previous work on the spectrum. To see this 
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vs. g = YMV /]vy/27r 3 / 2 for (a) AT = 5 (b) if = 6 and for various values of the transverse 
resolution T. 



behavior we study the large-r behavior of the correlator at fixed g as a function of the 
transverse resolution T and at fixed T as a function of the coupling g. We see a hint that 
something unusual is occurring in Fig. §. For values of the coupling up to about g = 1 
we see the typical rapid convergence in the transverse momentum cutoff; however, at larger 
coupling the convergence appears to deteriorate, and we see that for g = 5 the correlator is 
smaller than at g — 10. We see this same behavior at both K — 5 and K = 6.0 In Fig. § 
we see that the correlator does not in fact decrease monotonically with g but rather has a 
singularity at a particular value of the coupling which is a function of K and T. Beyond the 
singularity the correlator again appears to behave well. 

If we plot the 'critical' couplings, at which the correlator goes to zero, versus a/T, as in 

Fig. [5|, we see that they lie on a straight line, i.e. this coupling is a linear function of y/T 

in both cases, K — 5 and 6. Consequently, the 'critical' coupling goes to infinity in the 

transverse continuum limit. It appears as though we have encountered a finite transverse 

cutoff effect. The most likely conclusion is that our numerical calculation of the BPS wave 

function is only valid for g < g CT i t (T). While the large-r correlator does converge above the 

critical coupling, it is unclear at this time if it has any significance. It might have been 

expected that one would need larger and larger transverse resolution to probe the strong 

coupling region, the occurrence of the singular behavior that we see is a surprise, and we 

have no detailed explanation for it at this time. We see no evidence of a singular behavior 

2 We do not see this behavior at K = 4, but it is not unusual for effects to appear only at a large enough 
resolution in SDLCQ. 
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at small or intermediate r. This indicates, but does not prove, that our calculations of the 
massive bound states is valid at all g. 

We do not seem to see a region dominated the massive bound states, that is, a region 
where r is large enough that we see the structure of the bound states but small enough that 
the correlator is not dominated by the massless states of the theory. Such a region might 
give us other important information about this theory. 



6 Conclusion 



In this work we calculate the correlator (0|T ++ (x)T++(0)|0) in M = 1, SYM in 2+1 dimen- 
sion at large N c in the collinear limit. We find that the free-particle correlator behaves like 
1/r 6 , in agreement with results from conformal field theory. The contribution from an indi- 
vidual bound state is found to behave like 1/r 5 , and at small r such contributions conspire 
to reproduce the conformal field theory result 1/r 6 . We do not seem to find an intermediate 
region in r where the correlator behaves as 1/r 5 , reflecting the behavior of the individual 
massive bound states. 

At large r the correlator is dominated by the massless BPS states of the theory. We find 
that as a function of g the large-r correlator has a critical value of g where it abruptly drops 
to zero. We have investigated this singular behavior and find that at fixed longitudinal 
resolution the critical coupling grows linearly with \/T. We conjecture that this critical 
coupling signals the breakdown of SDLCQ at sufficiently strong coupling at fixed transverse 
resolution, T. While this might not be surprising in general, it is surprising that the behavior 
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appears in the BPS wave functions and that we see no sign of this behavior in the massive 
states. We find that above the critical coupling the correlator still converges but significantly 
slower. It is unclear at this time if we should attach any significance to the correlator in this 
region. 

This calculation emphasizes the importance of BPS wave functions which carry important 
coupling dependence, even though the mass eigenvalues are independent of the coupling. We 
will discuss the spectrum, the wavefunctions and associated properties of all the low energy 
bound states of Af = 1 SYM in 2+1 dimensions in a subsequent paper [19[ . 



A number of computational improvements have been implemented in our code to allow 
us to make these detailed calculations. The code now fully utilizes the three known discrete 
symmetries of the theory, namely supersymmetry, transverse parity P, Eq. (|15|), and the Z 2 
symmetry S, Eq. flUf). This reduces the dimension of the Hamiltonian matrix by a factor 
of 8. Other, more efficient storage techniques allow us to handle on the order of 2,000,000 
states in this calculation, which has been performed on a single processor Linux workstation. 
Our improved storage techniques should allow us to expand this calculation to include higher 
supersymmetries without a significant expansion of the code or computational power. We 
remain hopeful that porting to a parallel machine will allow us to tackle problems in full 
3+1 dimensions. 
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